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Abstract. The position operator (defined within Schrodinger representation as usual) 
becomes meaningless when the usual Born- von Karman periodic boundary conditions 
are adopted: this fact is at the root of the polarization problem. I show how to define 
the position expectation value by means of rather peculiar many-body (multiplicative) 
operator acting on the wavefunction of the extended system. This definition can be 
regarded as the generalization of a precursor work, apparently unrelated to the po- 
larization problem. For uncorrelated electrons, the present finding coincides with the 
so-called "single-point Berry phase" formula, which can hardly be regarded as the ap- 
proximation of a continuum integral, and is computationally very useful for disordered 
systems. Simulations which are based on this concept are being performed by several 
groups. 

INTRODUCTION 

In a phenomenological description of dielectric media the concept of macroscopic 
polarization is the basic one [1]. One formally defines the polarization of a macro- 
scopic sample as its electrical dipole moment, divided by the volume. The problem 
is that the dipole of a finite system is dominated by surface effects, which seem to 
make polarization ill defined as a property of the bulk. In condensed matter theory, 
one typically gets rid of undesirable surface effects by dealing with extended sys- 
tems within Born- von Karman periodic boundary conditions (BvK). Unfortunately, 
this does not solve the polarization problem: the dipole is in fact the expectation 
value of the quantum-mechanical electronic position, which in its usual form is an 
ill defined operator within BvK. 

Because of the above reasons, even defining what bulk dielectric polarization is 
remained a major challenge for many years. If one assumes the polarization el- 
ements as discrete, a la Clausius-Mossotti, then the polarization mechanism can 
be safely understood: but such oversimplified picture does not apply to a real di- 
electric, where the electronic distribution is continuous, and often very delocalized. 



Typically, textbooks attempt to explain polarization of a bulk periodic solid via the 
dipole moment of a unit cell, or something of the kind [2,3]. These definitions are 
incorrect [4]: according to the modern viewpoint, bulk macroscopic polarization is 
a physical observable completely independent from the periodic charge distribution 
of the polarized dielectric. 

The breakthrough was fostered by the 1992 Williamsburg meeting. Macroscopic 
polarization was defined in terms of the wavefunctions, not of the charge. This 
definition has an unambiguous thermodynamic limit, such that BvK and Bloch 
states can be used with no harm [5]. In the following months a modern theory of 
macroscopic polarization in crystalline dielectrics has been completely established 
[6,7], thanks to a major advance due to R.D. King-Smith and D. Vanderbilt [8], who 
expressed polarization in terms of a Berry phase [9,10]. An early comprehensive 
account of the modern theory exists [7]. Other less technical presentations are 
available as well [11,12]; for an oversimplified outline see Ref. [13]. 

In its original form, the Berry-phase theory of polarization was based on 
independent-electron wavefunctions, such as the Kohn-Sham orbitals of density- 
functional theory [7,14]. Its generalization to correlated many-electron wavefunc- 
tions is due to Ortiz and Martin [15]. Both the independent-electron and the 
correlated-electron versions of the Berry-phase theory of polarization rely on some 
lattice periodicity, and define macroscopic polarization in the form of a reciprocal- 
space integral. In particular, the Ortiz-Martin polarization of a correlated many- 
electron system takes the form of a peculiar "ensemble average", since it is an 
integral over a set of many-body ground states. This viewpoint is indeed cor- 
rect and very valuable: the Ortiz-Martin theory has even been implemented to 
study the polarization of interesting model systems of highly correlated electrons 
[16-18]. However, the fact that macroscopic polarization could not be even defined 
by means a "pure state" expectation value (over a single many-electron ground 
state) was a disturbing drawback. This has been overcome thanks to the advance 
of Ref. [19], where a novel solution to the polarization problem is provided. The 
compact formula arrived at in Ref. [19], Eq. (14) below, is apparently unrelated to 
the Berry-phase concept; it applies on the same footing to correlated systems and 
to independent-electron systems, as well as to crystalline and to disordered systems. 
Lattice periodicity and integration in reciprocal space are no longer needed in order 
to define what polarization is. Polarization can in fact be cast as the expectation 
value of a rather peculiar operator over the many-electron ground wavefunction (in 
the thermodynamic limit). 

In the present work, I present the main achievement of Ref. [19] under a rather 
different light than in the original paper. I will show that its main result can be 
regarded as the many-electron generalization of a precursor work, apparently unre- 
lated to the polarization problem [20] . I will also discuss in detail the independent- 
electron case, linking in particular the results of Ref. [19] to the so-called single- 
point Berry phase [21,22], a concept which proves very useful, particularly in the 
study of disordered systems. 



THE HILBERT SPACE OF CONDENSED MATTER 



Almost invariably, in condensed matter theory one adopts BvK for the electronic 
wavef unctions, when dealing with either crystalline systems or disordered systems. 
There are several good reasons for such choice. One reason is that one eliminates 
any surface by construction, thus getting rid of indesirable surface effects. Another 
reason is that in the crystalline case one can exploit the virtues of Bloch's theorem. 
I stress that the BvK choice is tantamount to defining the Hilbert space where our 
solutions of Schrodinger equation live. 

For the sake of simplicity, the present work will deal with the one-dimensional 
case. For a single-particle wavefunction BvK reads ifi(x + L) = ip(x), where L is the 
imposed periodicity, chosen to be large with respect to atomic dimensions. Notice 
that lattice periodicity is not assumed, and BvK applies to disordered systems 
as well. An operator maps any vector of the given Hilbert space into another 
vector belonging to the same space: the multiplicative position operator x is not a 
legitimate operator when BvK are adopted for the state vectors, since xvfj(x) is not 
a periodic function whenever ip(x) is such. Of course, any periodic function of x is 
a legitimate multiplicative operator: this is the case e.g. of the nuclear potential 
acting on the electrons. 

In the following, we will need to explicitly introduce many-body wavefunctions 
as well. BvK then imposes periodicity in each electronic variable separately: 

^ (x 1 ,...,x i ,...,x N ) = ^ (x 1 , Xi + L, ...,x N ). (1) 

Our interest is indeed in studying a bulk system: N electrons in a segment of length 
L. Eventually the thermodynamic limit is taken: L — > oo, N — > oo, and N/L = hq 
constant. We will also assume throughout the ground state nondegenerate, and we 
deal with insulating systems only: this means that the gap between the ground 
eigenvalue and the excited ones remains finite for L — > oo. 

In order to illustrate where the main problem is, let us start with an Hilbert 
space different from the BvK one, using instead the boundary conditions which 
are appropriate for the bound states of an isolated atom or molecule. In this 
case the many-body wavefunction $0(^1; • • • , x u ■ ■ ■ , x n) goes to zero exponentially 
whenever gets large. Within this Hilbert space there is no problem with the 
position, which is trivially defined as the multiplicative operator 

N 

* = E*i ( 2 ) 

1=1 

within the Schrodinger representation. The ground-state expectation value is then 
simply: 

(X) = ($ |X|$o) = J dxxn(x), (3) 

where n(x) is the one-particle density. The value of (X) scales with the system size, 
and the quantity of interest is indeed the dipole per unit length, which coincides 



with macroscopic polarization. The expectation value in Eq. (3) is dominated 
by surface effects, thus making the definition of a bulk quantity problematic. If 
instead one adopts BvK, Eq. (1), then the operator defined by Eq. (2) is no longer 
a legitimate one, as explained above. 

The major goal is therefore to define the expectation value of the electronic 
position (X) in the BvK Hilbert space, and to prove that our definition provides in 
the thermodynamic limit the physical macroscopic polarization of the sample: this 
has been very recently achieved in Ref. [19]. Before illustrating this major advance, 
we discuss an important precursor work, apparently unrelated with the polarization 
problem, where nonetheless the expectation value of the position operator plays the 
key role. 

THE ELECTRON-IN-BROTH FORMULA 

Some years ago, A. Selloni et al. [20] addressed the properties of electrons dis- 
solved in molten salts at high dilution, in a paper which at the time was commonly 
nicknamed the "electron in broth". The physical problem was studied by means 
of a mixed quantum-classical simulation, where a lone electron was adiabatically 
moving in a molten salt (the "broth") at finite temperature. The simulation cell 
contained 32 cations, 31 anions, and a single electron. KC1 was the original case 
study, which therefore addressed the liquid state analogue of an F center; other 
systems were studied afterwards [23]. The motion of the ions was assumed as com- 
pletely classical, and the Newton equation of motions were integrated by means 
of standard molecular dynamics (MD) techniques, though the ionic motion was 
coupled to the quantum degree of freedom of the electron. The electronic ground 
wavefunction was determined solving the time-dependent Schrodinger equation at 
each MD time step. Two snapshots of the simulation, reproduced from the original 
paper, are shown in Fig. 1. 

As usual in MD simulations, periodic boundary conditions are adopted for the 




FIGURE 1. Contour plot of the electronic density (integrated along the sight line) at two 
different time steps. The dots are projections of the ionic positions. After Ref. [20]. 



classical ionic motion. Ideally, the ionic motion occurs in a simulation cell which 
is surrounded by periodic replicas: inter-cell interactions are accounted for, thus 
avoiding surface effects. Analogously, the electronic wavefunction is chosen in the 
work of Selloni et al. to obey BvK over the simulation cell, and therefore features 
periodic replicas as well. In Fig. 1, the "center" of the electron distribution is close 
to one of the borders of the simulation cell: the contour plot clearly shows the tail 
of its replica spilling into the opposite border. 

One of the main properties investigated in Ref. [20] was the electronic diffu- 
sion, where the thermal ionic motion is the driving agent (within the adiabatic 
approximation). In order to perform this study, one has to identify first of all 
where the "center" of the electronic distribution is. Intuitively, the distributions in 
Fig. 1 appear to have a "center" , which however is defined only modulo the replica 
periodicity, and furthermore cannot be evaluated simply as / dr r\ip(r)\ 2 precisely 
because of BvK. Selloni et al. solved the problem by means of a very elegant and 
far-reaching formula. Here I write its one-dimensional analogue in the form: 

(x) = — Im log / dx e l2 ? x \ip(x)\ 2 . (4) 
2n Jo 

This formula provides the center of the electronic distribution, which is defined 
unambiguously though modulo L, as it must be. 



HOW THE FORMULA WORKS 

The density of the single-particle wavefunction n(x) = \ip(x)\ 2 obeys BvK peri- 
odicity and can therefore be written as a Fourier series: 



n 



(x)= £ C s e^». (5) 



The electron-in-broth formula, Eq. (4) becomes: 

L _ . r L . .-a*. , , L 



(X) = — 
X 1 2tt 



Im log / dx e l2 L x n(x) = — Im log C_i, (6) 
Jo 2ir 



a quantity which is well defined (modulo L) whenever the s — — 1 Fourier coefficient 
is nonvanishing. 

In general, any periodic function n can be written — in a nonunique way — as the 
superposition of a localized function n\ oc and of its periodic replicas. We write in 
the most general case the superposition as 



n 



( x ) = n ioc(x - x - nL), (7) 



where Xq is hitherto arbitrary. Then the Fourier coefficients of n(x) are related to 
the Fourier transform of n\ oc {x) as: 



We therefore express (x), Eq. (4), as: 



\x) = — Im log 



e^ Xo nio C (- — 



L _ 27T 

= 2:0 + — Im log ni oc (- — ). (9) 



We now specialize to the only interesting case: n\ oc (x) is a real positive function, 
which integrates to one over (—00,00); furthermore n\ oc (x) is centered at x = 
and decays over distances of the order of Ax. In the special case where n\ oc (x) is 
an even function, its Fourier transform is real and Eq. (9) simplifies to (x) = xq 
(mod L), which is indeed the expected result. In the general case where n\ oc (x) is 
not an even function, an expansion of its Fourier transform is needed. The Fourier 
transform n\ oc (k) decays over a reciprocal distance of the order of 1/Ax: therefore 
supposing Ax <C L (as e.g. is clearly the case in Fig. 1) the expansion yields: 

27j- roo 27T f°° 

n\oc( — -) = / dxn loc (x)+i— dxxn loc (x) + 0(L~ 2 ). (10) 

1j J —CO ±J J— OO 

Since by hypothesis ni oc (x) integrates to one, Eq (9) to leading order yields: 

/oo 
dxxn\ oc (x) (modL), (11) 
-00 

again a very meaningful expression for the center of the electronic distribution n(x). 



POLARIZATION AS A MANY-ELECTRON PROBLEM 

In the previous two Sections we have discussed how to define the expectation 
value of the position operator within BvK, when only a single electron is present 
in a large cell. When dealing instead with macroscopic polarization, we are faced 
with a system of many electrons, having a finite density. As anticipated above, our 
main interest is in studying a bulk system: N electrons in a segment of length L. 
Eventually the thermodynamic limit is taken: L — > 00, iV — > 00, and N/L = tiq 
constant. If one arrives at defining a suitable position expectation value (X), ana- 
logue of Eq. (3) for the much less trivial BvK case, then the electronic contribution 
to macroscopic polarization (dipole per unit length) is: 

Pa = lim (12) 

L^oo ±j 

where e is the electron charge. 

So the main issue is as follows: can the electron-in-broth formula, Eq. (4), be 
generalized to provide an expectation value {X) appropriate to the case of many 
electrons in the given BvK periodicity? The answer is yes, but a smart gener- 
alization is needed. The most straightforward generalization, where one simply 



identifies n(x) with the electron density of the extended system, does not work at 
all. There is a very good reason for this. 

Suppose the system is crystalline, and furthermore assume independent electrons 
for the sake of symplicity. According to the modern viewpoint [7] the macroscopic 
polarization is a bulk observable completely independent from the periodic charge 
distribution of the polarized dielectric. While the density is obtained from the 
square modulus of the Bloch orbitals, polarization can only be obtained from their 
phase. Indeed, if one attempts to actually evaluating Eq. (4) upon replacement 
of IV^)! 2 with the bulk density n(x) of the extended system a meaningless result 
is reached. In fact, if the lattice constant of the crystalline solid is a, the BvK 
periodicity is then taken over a multiple of a: L = Ma. Since n(x) is periodic with 
period a, its Fourier coefficients C s in Eq. (5) vanish unless s is a multiple of M: 
therefore C_i = 0, and the electron-in-broth formula, Eq. (6), is ill defined. 

In order to arrive at the correct many-electron generalization of Eq. (4), one first 
introduces the multiplicative many-body operator e l ~ x , where X is the simple 
operator defined in Eq. (2). In terms of this, one then defines: 



The similarity with Eq. (4) is self evident; as in the one-body case, the position 
expectation value is defined only modulo L. The main ingredient of our general- 
ized electron-in-broth formula is a many-body expectation value: an iV-electron 
integral where a genuine many-body operator appears. In general, one defines an 
operator to be one-body whenever it is the sum of N identical operators, acting on 
each electronic coordinate separately: for instance, the X operator itself is such. In 
order to express the expectation value of a one-body operator the full many-body 
wavefunction is not needed: knowledge of the reduced one-body reduced density 
matrix p is enough: this is e.g. the case in Eq. (3). I stress that, instead, the 
expectation value of e t ~ x over a correlated wavefunction cannot be expressed in 
terms of p, and knowledge of the iV-electron wavefunction is explicitly needed. The 
uncorrelated case where |^ ) is a single determinant will be discussed in the next 
Section. 

The electronic polarization, Eq. (12), becomes then: 



As anticipated in the Introduction, this compact and general expression for the 
macroscopic polarization applies on the same footing to correlated systems and to 
independent-electron systems, as well as to crystalline and to disordered systems. 
Notice that in any case L — > oo is an highly nontrivial limit, since the exponential 
operator in Eq. (14) goes formally to the identity, but the size of the system and the 
number of electrons in the wavefunction increase with L. Therefore some ad-hoc 
procedure is needed to demonstrate Eq. (14): this is reported in detail in Ref. [19], 
and will not be repeated here. 




(13) 




(14) 



INDEPENDENT PARTICLES AND SINGLE-POINT 

BERRY PHASE 



We focus here on an uncorrelated system of independent electrons, whose TV- 
electron wavefunction |^ ) is a Slater determinant. In this special case |^o) is 
uniquely determined by the one-body reduced density matrix p (which is the pro- 
jector over the set of the occupied spinorbitals): therefore the expectation value 
(X), Eq. (13), is uniquely determined by p. The expression for (X) — as well as 
the related one for the electronic polarization — is however most simply expressed 
directly in terms of the orbitals, as I am going to show. Such an expression, which 
goes under the name of "single-point Berry phase", was first proposed in a series 
of lecture notes [21]; it has been further scrutinized in some detail in a recent paper 
by Yaschenko et al. [22]. 

Suppose N is even, and |^o) is a singlet. The Slater determinant has thus the 
form: 

l^o) = ^=kl^l^ 2 • • • ¥N/2¥n/2\, (15) 

where fi are the single-particle orbitals. It is then expedient to define 

|* )= e ^|*o): (16) 

even |^o) is indeed a Slater determinant, where each orbital <fi(x) of |\&o) is mul- 
tiplied by the plane wave e l ~ x . According to a well known theorem, the over- 
lap amongst two determinants is equal to the determinant of the overlap matrix 
amongst the orbitals. We therefore define the matrix (of size N/2 x N/2): 

Sij= f L dx<p*(x)e i ^ x <p j (x) : (17) 

J 

in terms of which we easily get 

(X) = ^Im log (tfol^o) = -Im log det S, (18) 
2n 7r 

where the factor of 2 accounts for spin. 

Notice that nowhere we have assumed crystalline symmetry: the ipi orbitals will 
therefore have the most general form, while obviously obeying BvK at the sample 
boundary. Eq. (18) applies on the same ground to crystalline and to disordered 
systems: in the crystalline case it can be shown to be equivalent to a discretization 
of the usual continuum Berry phase [7], where crystalline Bloch orbitals are used. A 
simple test case is studied by Yaschenko et al. [22], where the numerical equivalence 
is actually demonstrated. 

The single-point Berry phase has been recently implemented quite successfully to 
calculate infrared spectra of disordered systems [24,25]. In a typical Car-Parrinello 



simulation [26], the thermodynamic limit is approximated using BvK over a large 
simulation cell, which contains several tens of atoms. The time evolution of the 
ions is followed along discrete time steps, while the electrons follow adiabatically. 
In order to evaluate the infrared spectrum, one needs to evaluate the electronic 
macroscopic current traversing the simulation cell. This is done on the fly with the 
use of the single-point Berry phase, where (X), Eq. (18), is evaluated at each time 
step, and its time derivative is approximated with a finite difference [24]. 

CONCLUSIONS 

The formula originally proposed by Selloni et al. [20] for dealing with a lone elec- 
tron in a large liquid sample, and named here the "electron-in-broth formula" , is 
quite naturally generalized to the many-body case. It provides then a novel solution 
to the polarization problem, Eq. (14): such solution, though apparently unrelated 
to the Berry-phase concept, is indeed a discretized Berry phase in disguise [19]. 

Since any Berry phase is by definition a phase, it may appear surprising that in 
Eq. (14) only the square modulus of the wavefunction enters: any phase information 
is then apparently obliterated. One has to bear in mind, however, that |^o) is the 
many-body wavefunction of the iV-electron system. In the independent-electron 
case, this wavefunction embeds the relevant information about the relative phases 
of the one-electron orbitals, and this is what one needs in order to evaluate the 
macroscopic polarization; in the correlated-electron case, the wavefunction carries 
in a more complex way the analogue information. 
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